Insulator-Metal Phase Diagram of the Optimally Doped Manganites 
from the Disordered Holstein-Double Exchange Model 
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We study the Holstein-Double Exchange model in three dimensions in the presence of substitutional 
disorder. Using a new Monte Carlo technique we establish the phase diagram of the clean model 
and then focus on the effect of varying electron-phonon coupling and disorder at fixed electron 
density. We demonstrate how extrinsic disorder controls the interplay of lattice polaron effects and 
spin fluctuations and leads to widely varying regimes in transport. Our results on the disorder 
dependence of the ferromagnetic T c and metal-insulator transitions bear direct comparison to data 
on the 'optimally doped', x = 0.3—0.4, manganites. We highlight disorder induced polaron formation 
as a key effect in these materials, organise a wide variety of data into a simple 'global phase diagram', 
and make several experimental predictions. 



The Holstein-Double Exchange (H-DE) model provides 
a minimal description of strongly coupled charge, spin, 
and lattice degrees of freedom, typical of the perovskite 
manganites Ai-^A^MnOs [1], and captures the inter- 
play of polaronic tendency and spin disorder that oc- 
curs in these materials. Although a detailed solution has 
not been available in three dimensions (3D), the H-DE 
model is believed to describe the trend towards stronger 
localisation and weakened ferromagnetism observed with 
reducing A site ionic radius (ta) in the manganites. 

The role of quenched disorder in these systems, aris- 
ing, for example, from the cation size mismatch (a a) is 
only partially understood. Most of the disorder related 
theory [2,3] has focused on effects near a first order phase 
boundary, i.e, the 'bicritical' region [4], motivated by the 
observation of cluster coexistence in some materials [5] . 
Phase competetion and bicriticality, however, occurs only 
in a limited part of the phase diagram. Studies with con- 
trolled variation in a a indicate that quenched disorder 
has a dramatic effect on the ferromagnetic T c and resis- 
tivity [6-9], and optical response [10] even far from bicrit- 
icality. The origin of these effects remain unexplained. 

Our primary contribution in this paper is to demon- 
strate how this unusual sensitivity arises due to a posi- 
tive feedback of quenched disorder on the polaron forma- 
tion tendency, and its interplay with spin fluctuations. 
Through a detailed real space solution of the disordered 
H-DE model on large 3D lattices we (i) demonstrate dis- 
order controlled metal- insulator transitions (MIT), in the 
ground state and at finite temperature, that are absent in 
the 'clean' system, (ii) explain the rapid suppression of 
T c with increasing disorder, (Hi) compare our results in 
detail with data on the manganites at x = 0.3 — 0.4, con- 
structing a 'global phase diagram', and (iv) make testable 
predictions about the a a dependence of optical spectral 
weight, lattice distortion, and tunneling density of states 
in the optimally doped manganites. 



We study the adiabatic disordered (d) H-DE model: 
H = -t>( c]„c,;„ + h.c) + 2j(e l - n)rii 
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The t are nearest neighbour hopping on a cubic lattice, 
ej is the quenched binary disorder, with e, = and value 
±A, A is the electron-phonon (EP) interaction, coupling 
electron density to the local distortion Xi, K is the stiff- 
ness, (U is the chemical potential, and Jh is the Hunds 
coupling. The parameters in the problem are X/t, A/t, 
electron density n, and temperature T. We treat the 
phonons and spins as classical (the adiabatic limit), as- 
sume Jn/t —> oo, and set |S| = 1. We also set iv = 1 
and measure energy, frequency, T, etc, in units of t. 

The clean H-DE model was proposed early on [11] as 
a minimal model for the manganites, and a mean field 
study [11] suggested strong decrease of T c and increasing 
localisation with increase in EP coupling. More recently 
the H-DE model with 'cooperative phonons' has been 
analysed [12] in 3D, yielding a thermally driven MIT 
at low n but not near optimality. The effect of disor- 
der has been considered on EP-DE models but with fo- 
cus mainly on phase coexistence [13] and bicriticality [3]. 
These studies yield insight on how disorder leads to clus- 
ter formation near a first order phase boundary [2,13], or 
its impact on a commensurate charge ordered phase [3], 
but do not have much bearing on the effects observed, 
for example, at x = 0.3, where there is no obvious phase 
competetion. In our understanding, the effects of disor- 
der in the manganites, away from bicriticality, await an 
explanation. In this paper we overcome a key technical 
limitation of the earlier studies and provide new results 
on transport combining the effects of strong EP coupling, 
thermal spin fluctuations, and quenched disorder, on lat- 
tices upto 12 3 in size. 



Since we work with Jjj/t — » oo, the electron spin 
is 'slaved' to the core spin orientation leading to an 
effectively spinless, hopping disordered model: H = 
-£<«>(*«7i7i + h.c)+^ i (e i -fi-Xx i )n i +{K/2)J2 i xl 
The hopping amplitude depends on the spin orientations 
via Uj/t = cosfjcos-^- + sin^sin^-e 1 ^^^^, where 9i 
and <f>i are respectively the polar and azimuthal angles 
of the spin Si. We use our recently developed "travel- 
ling cluster approximation" (TCA) [14] for solving the 
annealing problem on large lattices, by computing the 
energy cost of a MC update by constructing a 'cluster 
Hamiltonian', of N c — L\ sites, around the update site. 
We will use iV c = 4 3 and, as documented earlier [14], this 
is sufficient for equilibriating the spin and lattice vari- 
ables. The electronic properties, e.g, resistivity, density 
of states, and optical response, are obtained by diago- 
nalising the full electron Hamiltonian on the equilibrium 
{xi,Si} configurations, averaging thermally and finally 
over 8 — 10 realisations of disorder {ei}. The transport 
results are obtained using a formulation described in de- 
tail previously [15]. We will use system size N — 8 3 , but 
have checked our results on N upto 12 3 . Our resistivity 
results are in units of hag /(ire 2 ), where ag is the lattice 
spacing. As a rough estimate, p = 1 in our units would 
correspond to ~ 30/iilcm. 

Let us summarise the phases of the clean H-DE model 
before addressing the effects of disorder, (a). Struc- 
turally, a phase can be either charge ordered (CO), or 
charge disordered, i.e a polaron liquid (PL), or homoge- 
neous, i.e & Fermi liquid (FL). (b). Magnetically it can 
be either ferromagnetic (F) or paramagnetic (P), and 
(c). in terms of electrical response it can be a 'metal', 
with dp/dT > 0, or an insulator, with dp/dT < 0. 

Combining these features, the ground state can be ei- 
ther a ferromagnetic Fermi liquid (F-FL), or F-PL, or 
F-CO, depending on A and n, while finite T allows P-FL, 
P-PL, and P-CO phases as well. Notice that while the 
threshold for single polaron trapping, i.e n = 0, in 3D, is 
E p = \ 2 /2K « 5.4i, i.e X c (n = 0) ~ 3.3 [16], at finite 
density the threshold for a collective PL is lower, e.g, 
X c (n = 0.3) ~ 2.6. These thresholds are further affected 
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FIG. 1. Phase diagram in the clean problem: (a) and (b) 
show the phase diagram at A = 2.0 and A = 2.5 respectively. 
The notation for the phases is explained in the text. The 
chessboard pattern is for a charge ordered (CO) phase with 
(7r,7r,7r) modulation. Shaded regions represent coexistence. 



by structural and spin disorder. 

For A < 2.0, Fig.l.(a), the phases arc (i) F-FL at low 
T away from n = 0.5, (ii) P-FL at high T, and (Hi) a 
F-CO phase at n = 0.5, separated from the F-FL, at low 
T, by a window of phase separation. For A = 2.0 the 
T c at n = 0.5 is suppressed with respect to the pure DE 
problem, but is still larger than the CO scale Tqo- At 
stronger coupling, of which A = 2.5 is typical, Fig. 1.(6)., 
much of the phase diagram is taken up by insulating 
phases. As A is increased beyond 2.5 the metallic win- 
dow at low n shrinks quickly and all the phases turn 
insulating. The ground state, however, continues to be 
ferromagnetic since we have not included any competing 
antifcrromagnctism (AF). From the resistivity we con- 
firmed that the clean model does not show a thermally 
driven MIT except for n < 0.1, as noted earlier [12]. 

Even weak structural disorder dramatically affects the 
electronic state at strong coupling, leading to metal- 
insulator transitions which are not present in the clean 
limit. In the rest of the paper we focus on n = 0.3, typ- 
ical of the high density metallic regime, and study the 
effect of varying disorder and EP coupling. Since the ta 
controlled bandwidth variation in real materials is rela- 
tively small we use only a ±10% variation in X/t about 
A = 2.0 and study the weak disorder regime, A = — 1.0. 

Disorder creates density inhomogencities in the FL 
which are amplified by the EP coupling, leading to a 
rapid rise in residual resistivity, and also opens up a win- 
dow for the crucial F-FL to P-PL crossover. Fig. 2. (a) — 
(c) show p(T) with varying A and A illustrating the three 
broad transport regimes. At A = 1.8, which is towards 
the lower end of our coupling range, moderate A in- 
creases p(0) but the overall p(T) has the same nature 
as in the simple DE model [15] and neither EP coupling 
nor disorder seems to have any qualitative effect. This 
is true of A < 0.6 and, we suggest, is akin to what 
one observes in Lao.7Sro.3Mn03. At A = 1.0, where 
p(0) > 100 ~ pMott w e see hints of a weak downturn 
in p(T) at high T. Overall, p(T) at A = 1.8 is mainly 
characterised by a F-FL to P-FL crossover. The density 
of states (DOS), Fig.2.(d), remains featureless at all T 
and the lattice distortions, Fig.2.(g), are weak with no 
significant variation across T c . This "DE like" charac- 
ter in p(T) survives upto A = 2.2 in the clean limit, see 
Fig. 2. (6) — (c). However, for A > 2.0 addition of weak 
disorder (i) sharply increases p(0) and creates a pseudo- 
gap in the DOS, (ii) maintains a regime with dp/dT > 
for T < T c , but (Hi) crosses over to an 'insulating' PL 
regime with dp/dT < for T > T c . As Fig.2.(/i) - (i) 
indicate, these changes are associated with large lattice 
distortions and their variation across T c . For A ~ f.O 
the ground state itself turns insulating for A > 2.0. We 
provide a compact picture of the various MIT in Fig. 3, 
discussed later. 

To gain insight on the F-FL to P-PL crossover we also 
calculated the 'effective carrier number', n e ff(u),T) = 



L a(u>, T)duj, at T = and ui — 0.5, as well the density 
field n(r. T). The combination of n e //, Fig. 4. (a), the spa- 
tial patterns (not shown here) and the DOS, Fig.2.(e), 
suggest the following scenario at strong coupling and 
moderate disorder: (i) there is polaronic localisation of 
a fraction of carriers (states at u> < p) at T = 0, lower- 
ing n e ff , but the states at the chemical potential remain 
extended, (ii) the polaronic states are further localised 
with increasing T due to DE spin disorder, while the de- 
localised electrons near p experience scattering from the 
spin fluctuations, leading to dp/dT > 0, (Hi) by the time 
T ~ T c the net disorder arising from ej, the lattice dis- 
tortions, and spin disorder is so large that the system 
is effectively at a localisation transition on the complex 
{ei, Xi, Sj} background: the states at p also get localised, 
(iv) the T > T c phase has a modest 'mobility gap' (at 
moderate A, A) so there is no simple activated behaviour. 
The pseudogap in the DOS, Fig.2.(e), deepens initially 
and gradually fills up for T > T c . 

Let us consider the relevance of these results to man- 
ganites at s ~ 0.3 — 0.4, away from commensurate fill- 
ing and charge ordering instabilities. The mean hopping 
amplitude between Mn ions, via the oxygen orbitals, is 
controlled by ta through the Mn-O-Mn bond angle. The 
mismatch between different A site cations is quantified 
through the variance, a a, of the ionic size, and acts as 
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FIG. 2. Colour online: Resistivity, p(T), DOS, and lat- 
tice distortion at n = 0.3. (a) — (c). p(T), with varying A 
and A. (d) - (/). DOS, N(u,T) for the indicated parame- 
ters, (g) — (i). RMS lattice distortion, 5x rrns — \/{(Xi — x) 2 ), 
where x is the system averaged distortion. 
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FIG. 3. Colour online: Transport regimes and ferromag- 
netic T c : comparing data on the x = 0.3 manganites with 
our n = 0.3 results, (a). Transport 'phase diagram' indicat- 
ing the nature of low T to high T crossover with changing ta 
and (ta, data from Ref[6-9] (6). T c (ga, ta) in the manganites, 
from Ref[6-7], (c) transport crossovers from our calculation, 
and (d). our results on T C (A, A). 

disorder on the electronic system. For a material of com- 
position Ai-^A^MnOa, say, ta and a a depends on the 
radius of the A, A' ions as well as x. To disentangle the 
effect of varying carrier density (x) from that of varying 
'coupling constants', there have been systematic studies 
of the manganites at fixed x, with controlled variation in 
ta, related to our X/t ratio, and a a, related to our A/t. 

Focusing initially om = 0.3, the data from a wide vari- 
ety of sources can be organised in terms of the transport 
response, p(T), and magnetic T c . We choose a simple 
characterisation where we construct the MIT phase dia- 
gram in terms of the low T to high T crossover observed 
in p(T). For example, Lao.7Sro.3Mn03 (LaSr), is a metal 
at both low T and T > T c so it shows M-M crossover. In 
that spirit LaCa shows M-I crossover, while PrCa is I-I. 
These simple 'binary' cationic systems are well studied, 
have relatively small a a, and the transport response can 
be argued to arise mainly from variations in r a- However, 
more complex cationic combinations have been used, no- 
tably by Attfield and coworkers [6,7], to vary a a at fixed 
ta and x. Our organisation of the experimental data, 
Fig. 3. (a) uses three such families, at ta = 1-26, 1.23 and 
1.20A, as well as data from other sources [8,9] to highlight 
the combined impact of ta and a a variation on transport 
character. We have used the 9 fold coordinated values of 
ionic radius from Shannon [17] in computing r a and a a 
for about 30 compositions. 

Four features are notable in the experimental A4IT and 
T c data, Fig. 3. (a) — (6), at x = 0.3: (i) the 'clean' limit 
allows only a narrow window, ta ~ 1.18 — 1.20A, for M-I 
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FIG. 4. (a) Computed optical spectral weight, n e ff, at cut- 
off u = 0.5£ and T = 0.01 for varying A and A. Inset: data 
from Ref[10] on LaSr, etc, at x — 0.4, to — O.leV, plotted ver- 
sus oa- The ta values, LaSr downwards, are 1.256 A, 1.232 
A, 1.22 A and 1.202 A. (b) Our results on disorder induced 
lattice distortion at T — 0.01. Inset shows the a a induced in- 
crease in the mean square oxygen displacement at low T in a 
sequence of x = 0.3 and ta = 1.23A manganites, from Ref[6]. 

crossover, (ii) the M-I region broadens with increasing a a 
at the cost of the M-M region (red/dot-dash line), disor- 
der can induce a thermally driven MIT, {Hi) at T = just 
increasing a a can drive a metallic ground state insulating 
(blue/dashed line), (iv) the same a a has much stronger 
impact on transport and T c at small ta compared to 
large ta, for example a a ~ 0.01 A 2 has only modest ef- 
fect on T c at rA — 1.26 A, but leads to a large reduction 
at rA i$ 1.23 A. While we cannot make a direct numerical 
correspondence between a a and A/i, etc, our results on 
transport crossovers, Fig.3.(c), and T c , Fig.3.(d), show 
that all the trends (i) — (iv) above are correctly repro- 
duced. A fit to the detailed T dependence of resistivity 
will of course require a more sophisticated phonon model. 
Our T™ ax ~ 0.12* - 300K (for * - 0.25eV) compares 
reasonably with the maximum in Fig. 3. (b). 

Although the similarity of our results to experimental 
data on the manganites seems persuasive, let us exam- 
ine how these results may be modified in a more realistic 
model, (i) The manganites involve two e g orbitals cou- 
pled to John-Teller (JT) phonons, with the possibility of 
orbital order. Our preliminary calculations with the 2 
band JT model in 2D [18] suggests that for x ~ 0.3 - 0.4 
there is no non trivial orbital order and the transport re- 
sponse is similar to what we observe here. The one band 
model does not capture the correct physics as x — > 0, 
or x = 0.5, but appears reasonable for x ~ 0.3 — 0.4. 
(ii) The oxygen atoms are corner shared among Mn06 
octahedra so the phonon variables are cooperative. We 
have not explored this aspect but believe the main effects 
would be lower critical EP coupling for localisation [12] 
and possibly a sharper MIT across T c . (Hi) There are AF 
interactions present in the manganites. We repeated the 
d-H-DF, calculation in the presence of AF coupling and 
discovered that it leads to a sharper drop in T c , shifts the 
MIT phase boundaries somewhat, and leads to an AF-I 
phase at large A or A. The qualitative results of this 



paper, however, remain unaffected. 

Finally, if our argument about disorder induced po- 
laron formation is correct then the following features 
should be experimentally observed with increasing a a'- 
(a) systematic suppression of low frequency optical spec- 
tral weight, n e fj at u) < O.leV, at low T, (b) increase in 
the incoherent lattice distortion (oxygen Debye- Waller 
factor), progressively larger at smaller rA, and (c) ap- 
pearance of an increasingly deeper pseudogap in the tun- 
neling DOS at low T. Fig. 2 and Fig. 4 quantify these 
effects. Experiments in the recent past do claim disorder 
induced polaron formation in Lao.54Bao.46Mn03 [19] but 
more direct confirmation for x ~ 0.3 is desirable. 

To conclude, wc have studied a model with Holstcin 
phonons, double exchange and quenched disorder, and il- 
lustrated the disorder enhancement of polaronic tendency 
and its interplay with spin fluctuations. Our results com- 
pare well with the disorder dependence of transport and 
T c observed in the 'optimally doped' manganites. Wc 
make three testable predictions regarding the behaviour 
of optical spectral weight, oxygen Debye- Waller factor 
and tunneling density of states for controlled variations 
of a a in these materials. 

We thank HRI for use of the Beowulf cluster and 
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